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Classical light bending is investigated for weak gravitational fields in the presence of hypothetical 
local Lorentz violation. Using an effective field theory framework that describes general deviations 
from local Lorentz invariance, we derive a modified deflection angle for light passing near a massive 
body. The results include anisotropic effects not present for spherical sources in General Relativity 
as well as Weak Equivalence Principle violation. We develop an expression for the relative deflection 
of two distant stars that can be used to analyze data in past and future solar-system observations. 
The measurement sensitivities of such tests to coefficients for Lorentz violation are discussed. 
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I. INTRODUCTION 

A classic prediction of General Relativity (GR) is the 
bending of distant starlight by the Sun [1]. This was 
first confirmed in 1919 by Dyson, Eddington, and David- 
son to agree with Einstein's calculations of 1.75" at one 
Solar radii to within 30% accuracy [2]. More recent op- 
tical measurements during solar eclipses have made only 
marginal improvements [3]. The inclusion of radio as- 
tronomy has increased the accuracy of light deflection 
measurements to within 0.02%, providing additional firm 
evidence for the validity of the light deflection predicted 
in GR [4]. Measurements of the closely-related Shapiro 
time delay have also seen vast improvements recently. 
The analysis of two-way radio tracking of the Cassini 
probe matched the predictions of GR to within parts in 
100,000 [5]. 

Although it is currently the best fundamental theory 
of gravity, there remains widespread interest in develop- 
ing more precise tests of GR, including improved mea- 
surements of the bending of light, among others. These 
efforts are in part motivated by the intriguing possibility 
of finding deviations from GR. Such deviations could be a 
signature of a more fundamental unified theory of physics 
that successfully meshes GR with quantum theory and 
the Standard Model of particle physics. 

One possible signature that has been sought in many 
sensitive tests are minuscule violations of local Lorentz 
invariance, a fundamental tenet of GR [6]. Theoretical 
scenarios in which local Lorentz symmetry could be bro- 
ken are currently numerous in the literature, with early 
motivation coming from string field theory [7]. 

In order to investigate violations of local Lorentz in- 
variance, it is useful to have a theoretical framework in 
which to report measurements. One systematic frame- 
work for studying signals of Lorentz violation employs 
effective field theory. The idea is to incorporate known 
physics from GR and the Standard Model of parti- 
cle physics, into an effective action that also includes 
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generic Lorentz-violating terms. The additional Lorentz- 
violating terms in the action are controlled by coefficients 
for Lorentz violation, which are general coordinate ten- 
sor quantities describing the degree of Lorentz violation 
for each type of interaction (gravity, electrodynamics, 
etc.). These coefficients can be thought of as effectively 
fixed background fields in spacetime that couple to cur- 
vature and matter fields, though their origin can be dy- 
namical [7, 8]. The framework constructed in this man- 
ner is known as the Standard-Model Extension (SME) 
[9, 10], and has been adopted for numerous tests involv- 
ing light, matter, and gravity [11]. Connections between 
this framework and various classic test models for Special 
Relativity are discussed in Refs. [12]. 

Our focus in this work is on the signatures of Lorentz 
violation for gravitational tests. In the gravity sector 
of the SME, key signals in a number of experiments 
and observations have been established in Refs. [13-17]. 
Measurements constraining the coefficients for the grav- 
ity sector have already begun using atom-interferometric 
gravimetry [18, 19], lunar laser ranging [20], and short- 
range gravity tests [21]. In this paper, we analyze one 
of the fundamental tests of GR, the bending of light, 
in the effective field theory framework of the SME. This 
complements recent work on the related time-delay effect 
[14, 15]. 

We begin by deriving a general formula for the de- 
flection angle in section II A in terms of an arbitrary 
post-newtonian metric. The post-newtonian metric is 
described in Sec. II B. Assuming a stationary point-like 
mass we obtain the deflection angle in a limiting case in 
Sec. IIC, and a more accurate expression in Sec. II D. In 
Section III, we apply these results to light-bending tests 
in the solar system. We develop an expression for the 
measurable angle between two stars in Sec. Ill A. Details 
of the relative deflection angle and methods of analysis 
are discussed in Sec. IIIB. We illustrate the observable 
signals for Lorentz violation using a near-conjunction ex- 
ample in Sec. IIIC. Finally, in Sec. IV, we summarize 
the work and estimate the potential measurement sensi- 
tivities of existing and future light-bending tests to the 
coefficients for Lorentz violation in the gravity sector. 
Throughout this work we adopt notation and conventions 
as contained in Refs. [13 15]. In particular, we work in 
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natural units where c 
metric has signature - 



1 and the Minkowski spacetime 



II. THEORY 



A. Deflection basics 



The deflection a is the shift in the direction that light 
propagates from a straight line trajectory. We adopt a 
simplified gravitational lensing or light-bending scenario 
that involves a source S, a mass called the lens L, and an 
observer O. The geometric optics limit of electrodynam- 
ics in curved spacetime is assumed [24]. For a point-like 
lens the apparent source position observed by O is S a 
(see Fig. 1) [22, 23]. We assume the lens L, the source S, 
and the observer O are stationary throughout the light 
ray's propagation. The light ray emission is the event 
with coordinates (t e ,r J e ) and the observation of the light 
ray has coordinates (t p ,r p ). 




where dx 3 /dl is the tangent vector to the light path. The 
unit vector h 3 is the direction of the zeroth-order tangent 
to the light path. The light trajectory spatial endpoints 
are r 3 e and r 3 , which correspond to the I parameter values 
l p and l e . Referring to Fig. 1, l e = —r e ■ h and l p — r p - n. 
It is useful for later calculations to complete the set n 
and b with a perpendicular unit vector called f defined 
by f = h x b. 

If we apply the variational form of Fermat's principle 
(1) using the effective index of refraction (3) we obtain 
equations of motion for the light ray: 



dl 2 



5 jk 



dx 3 dx k 
~dl~dl 
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dx l 
~dl 



(4) 



Note that the terms on the right-hand side are perpendic- 
ular to dx 3 /dl, consistent with the definition of euclidean 
arclength I. This equation is equivalent to the geodesic 
equation for light to post-newtonian order v 2 , or PNO(2). 

The deflection a 3 follows by integration of Eq. (4) from 
the distant source S, at position r J e and I = —l e , to the 
observer O at position r p and I — l p . To PNO(2), the 
resulting deflection is given by the expression 



h 0j + 



dx 
~dJ 



>-kj 



(5) 



FIG. 1: The basic light bending scenario. The vector f p ex- 
tends from the central body L to the observer, f e is the vector 
from the central body to the emitter, n is in the direction of 
the unperturbed path, and b is the impact parameter vector 
(perpendicular to n). 

To calculate the deflection we can exploit Fermat's 
principle: the null geodesic path from (t e ,r{) to the ob- 
server's worldline is equivalent to the extremization of 
the arrival time t on the observer's worldline [25]. For 
a stationary observer, Fermat's principle is equivalent to 
the variational principle: 

5 jndi = 0. (1) 

Here n is the effective index of refraction of the gravi- 
tational field, I is the euclidean arclength {dl = V dx 2 ), 
and they are related by dt = ndl. 

The spacetime metric is expanded around a 

Minkowski background r/^ according to 

g^ v = Vn» + V- ( 2 ) 

Using the null condition for a light ray (g fJlU dx fl dx l/ = 0) 
we can evaluate n to leading order in metric fluctuations 

as 

i t . dx 3 -i , dx 3 dx k . 
n^l + gfooo + ft 0j -— + gfya , (3) 



where the symbol _L indicates a projection perpendicu- 
lar to dx 3 /dl, as in Eq. (4). The first integral in (5) 
is evaluated using the Euclidean arclength dl along the 
zeroth-order direction of the light ray. The second term 
is to be evaluated at the endpoints and plays a role in the 
result for the case where the observer is near the massive 
body L. 

The result in Eq. (5) applies to any metric that can 
be expanded around a Minkowski background in a post- 
newtonian series, so long as light behaves conventionally 
in the geometric optics limit (i.e., light follows a null 
geodesic). In the sections that follow, we shall apply 
this result to the post-newtonian metric that incorpo- 
rates local Lorentz and WEP violations using the SME 
framework. 



B. Post-newtonian metric 

We focus on the dominant terms in the gravitational 
sector of the SME framework [10]. This includes terms 
augmenting the pure-gravity sector (terms amending 
the Einstein-Hilbert action) as well as Lorentz-violating 
terms arising from the matter action. For the case of 
linearized gravity, the leading corrections to the post- 
newtonian metric of GR have been established and are 
discussed in detail in Refs. [13, 15]. Using a convenient 
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choice of coordinates and existing constraints on the vac- 
uum birefringence of light [26] , we can ignore Lorentz vi- 
olation in the electromagnetic sector, and hence assume 
light propagates normally [15]. 

The Lorentz violation for the pure-gravity sector is 
controlled by 9 coefficients called ~s^ v . For the matter sec- 
tor, the relevant coefficients are (a e s )n and c^„. The com- 
bined PNO(2) metric in harmonic coordinates is given by 



.9oo 

9jk 



-1 + [2 + 3^00 + 2cf + 4^(af ff ) ] U + s jk U jk , 
fry + %<j&)jP + isok + %(a s ctf ) k ]W k , 
5 jk + [2 - s 00 + 2qg, - 2%{a%B)o]SjkU 

+SjkSl m U 1 ' 

+ [2s 00 + 2^(a s eS ) }W k . 



SjiU lk -s ki U lj 



(6) 



/fi and Cqo, de- 



The coefficients for the matter sector, (a^ s ) 
pend on the structure of the source body S with mass to, 
and therefore will differ for distinct source bodies. Thus 
(af ff )^ and Cq indicate the presence of apparent Weak 
Equivalence Principle (WEP) violation as well as Lorentz 
violation. Note, however, that these WEP-violating coef- 
ficients do not affect the propagation of light rays, except 
through the spacetime metric. Also, as discussed in Ref. 
[15], a model dependent scaling a appears multiplying 
the coefficients (af ff ) M . In contrast, the coefficients s Mi , 
do not depend on the nature of the source body and 
therefore describe WEP-conserving local Lorentz viola- 
tion for gravity. 

Some limiting cases of this metric should be noted. 
This post-newtonian metric can be viewed as enlarging 
and complementing the Parametrized Post-Newtonian 
(PPN) metric [27] as discussed, for the case of vanish- 
ing matter coefficients, in Ref. [13]. Also, in the limit 
that the coefficients s M „, (of ff ) M , and Cq vanish, this re- 
sult reduces to the post-Newtonian metric of GR. 

For classical light bending scenarios we assume the 
source body L to be a point-like mass M. Furthermore, 
the source body is taken as the origin of the coordinate 
system. Thus, the dominant contributions to the poten- 
tials U and f7 jfe depend on the position of the test body 



U 

U° k 



GM 
r 

GMr^r k 



(7) 
(8) 



where G is Newton's gravitational constant. 



C. Light deflection: grazing case 

As a preliminary investigation of the modifications to 
the standard GR light-bending formula, we consider the 
simplified "grazing" case. In this scenario both the emit- 
ter and receiver are assumed to be very far from the 
source, thus we take l p — > oo and l e — > oo. Using the 



metric (6) in the general result (5) we obtain for this 
limiting case the deflection 



-4GM 



[l + *oo + c& + #(<&)<> 



+ (sofc + f (a S cS ) k )n k ]V - s kl b k T l ri 



(9) 



The standard GR result is obtained in the limit that the 
coefficients for Lorentz violation vanish. 

To illustrate some of the features of the result (9) we 
employ vector field plots indicating the apparent shift 
of incoming starlight as the source appears in front of a 
background of initially uniformly distributed stars. This 
is similar to methods that have been adopted for the GR 
signal in Refs. [28, 29]. For this simulation we focus on a 
small patch of sky centered on the deflecting body. We 
use x and y coordinates on the patch to make a vector 
field plot of different parts of the result (9). These co- 
ordinates will lie in the plane perpendicular to the unit 
vector n, which we can take to point (approximately) 
from the deflecting body to the observer. 1 

The left panel of Fig. 2 shows the initial star field in 
the absence of the source. For simplicity, we assume 
the distribution of stars is uniform. The right panel 
of Fig. 2 shows the deflection in the GR limit. In the 
expression for the deflection (9), the term proportional 
to V is scaled by the coefficients for Lorentz violation 
soo + Coo + ( a /-^0("off)o an d the n-dependent combi- 
nation (s fe + ( a /M)(a^ s )k)n k . Note that, due to the 
appearance of these coefficients, this scaling will change 
depending on the location of the patch of sky considered 
and the internal nature of the source. That this occurs 
is due to the anisotropy of the coefficients for Lorentz 
violation and their composition dependence. 

The last term in Eq. (9) points in the direction f , or- 
thogonal to b. This term is proportional to the combina- 
tion of coefficients s k ib k r l . If we express this combination 
in terms of the (x, y) coordinates on the patch of sky con- 
sidered, this can be expanded into the two independent 
combinations s xx — s yy and s xy . The deflections due to 
these two combinations of coefficients arc plotted in Fig. 
3, where we have set their values to one for illustration. 
These plots indicate that anisotropic deflection occurs for 
nonvanishing coefficients Sj k . In particular, this means 
that a light ray would deflect off of an otherwise flat 
plane defined by n and b, contrary to the conventional 
GR point-mass case. Note that these anisotropic effects 
can in principle be distinguished from the GR deflection 
due to higher multipole moments of the deflecting body 
[29, 30] by virtue of their 1/6 dependence. 



1 Strictly speaking, n points from the position of each star to the 
observer, but in the "grazing" approximation we can ignore the 
variation over the patch of sky. 
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FIG. 2: Initial uniform star field (left). Apparent shift of 
distant light in the standard general relativity case (right). 
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FIG. 3: Anisotropic apparent shift of star field due to the 
s X x — s yy coefficients (left) and the s xy coefficients (right). 
The local x coordinate runs horizontally and y is vertical. 



It is interesting to compare this result with that ob- 
tained in Rcfs. [14, 15] for the gravitational time delay 
in the presence of local Lorentz violation. The vector co- 
efficients soj and (a/M)(a^ s )j appear in the light bend- 
ing result (9), while they are absent in the round-trip 
time delay signal. Also, light-bending observations that 
involve one orientation of the observer, deflecting body, 
and the distant starlight could gain access to sets of coef- 
ficients distinct from those in a dedicated time-delay test 
with a different orientation of the relevant bodies. Note 
that this feature does not occur in the PPN formalism, 
for which both types of tests access the same parameter 
regardless of the underlying orientation. 



D. Light deflection: general case 

For applications of (5) to observations, the light- 
grazing approximation previously assumed must be re- 
considered. Here we still assume l e — > 00, but l p is 
treated as a relevant, finite term. This corresponds to 
the case where the observer is a finite distance from the 
lens L, while the light source is effectively at spatial in- 
finity. Referring to Fig. 1 it will be useful to define an 
angle $ between f p and n (thus r p ■ n = cos<I>). 

Evaluating the deflection formula (5) using the metric 



(6) for the case l e — > 00, we obtain the deflection 
-GM 



+ (5 0fc + ^(c4) fc )n fe )(l + cos$) 

+ (^00 + iHafff )o - s kl n k h l ) cos $ sin 2 $ 

+ 2(s ofc + § (a^ s )k)b k sin$ - s k in k b l sin 3 $j 
[2(s ok + ^(a s cS ) k )T k sin $ + s kl n k T l sin 3 $ 



- s ki b k r (2 + 2 cos $ + cos $ sin 2 $) 



(10) 



where the substitutions of l p /r p — cos $ and b/r p = sin $ 
have been made. Note that in the light-grazing limit the 
sin $ term vanishes and the cos $ term approaches unity, 
which results in the previous deflection angle presented 
in (9). 

The result (10) indicates that additional projections of 
coefficients for Lorentz violation arise in the more accu- 
rate deflection formula. Also, these additional combina- 
tions of coefficients appear to be distinct from the those 
that occur in the gravitational time delay derived in Refs. 
[14, 15]. By itself, the deflection in Eq. (10) is not directly 
measurable. This is because only the apparent position 
of a given source star is actually measured (at least dur- 
ing a single observation period). Instead, a comparative 
measurement is needed based on two or more observa- 
tions. In the following section we apply this result to 
calculate the relative deflection of two stars, which is a 
measurable quantity. 



III. SOLAR SYSTEM TESTS 

The key observable of interest for typical light bending 
measurements is the relative deflection of two (or more) 
stars. We seek here the angle between a "source" star 
and a "reference" star. In what follows quantities asso- 
ciated with the reference star will have r subscripts, and 
those for the source star will have no subscripts. The 
apparent positions of both of these stars will in principle 
be gravitationally deflected if the lens is near the line of 
sight. By continuous monitoring of the relative position 
of these two stars, one can measure the effects of the de- 
flection (10). Furthermore, the formula we derive below 
can in principle be applied repeatedly to systems of many 
stars. 



A. Relative deflection 

To calculate the relative deflection we adopt standard 
methods in the literature [27, 31]. We begin with a gen- 
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eral coordinate invariant expression for the angle VP be- 
tween two stars: 



cos * = 1 



P^Pr) 



(11) 



In this expression p^ and (p r ) M are the tangent four- 
vectors for the source light ray and the reference light 
ray, respectively. The four- velocity of the observer mea- 
suring the relative angle is u M . Evaluating this expression 
to PNO(2), and neglecting aberration terms, we obtain 



cos ^ = h ■ h r + h 
where a e g is given by 



+ n r ■ a cf f, 



"eft = a 3 + a(" '' h kj)±]t=i p 



(12) 



(13) 



The deflection a 3 is obtained from Eq. (10) and quantities 
with the label r are obtained by re-labeling all quantities 
occurring in the expression, involving the direction of the 
light ray, with the subscript r (e.g., h — > h r , b — > b r , etc. 
). The zeroth order, or straight line, trajectories from 
the source and reference stars are depicted in Fig. 4. 




FIG. 4: Solar light-bending scenario in the presence of a ref- 
erence source. For the source and reference star, we have the 
impact parameters 6 and b r , and directions h and n r , respec- 
tively. The reference plane shown is that spanned by b and 



We define an angle that represents the unperturbed 
angle between the two stars: 



n = cos 



(14) 



Thus is the angle between the two stars in flat space- 
time in the absence of other conventional effects such as 
aberration. The shift or change in the angle between the 
two stars is defined to be 



Since we are working in the post-newtonian approxi- 
mation it suffices to assume 5^ is small so that cos ^> ss 
cos — 3^ sin Using this approximation, the previ- 
ous result (10), and the metric in (6), we obtain 

sin * 5* ~ (h r ■ b[2(l + cos $) + B] + n r ■ tT) 

+ ~T~~ (" ' ^'t 2 *- 1 + C0S< M + B r ] + h ■ T r T,^j 

(16) 

In (16) the terms B, T, B r and T r are proportional to 
combinations of the coefficients for Lorentz violation s 



and (af ff ) M . Explicitly, they are given by 



B 



2[s 00 



( -oo 



fH a fff)o 



B r = 



T 



+ (sofc + f (a c S ff )fe)n' c ](l+cos$) 
+ \{s k ib k b l - Skih k h l ) cos $ sin 2 $ 
+2(s fe + §(a^ s ) k )b k sin$ 

- Skih k b l sin 3>(sin 2 $ - ±) 

2[s o+4> + # (a c 5 ff ) 

+ (s fc + #(af ff )fc)n£](l + cos$ r ) 
+ \{'skib k b r - s k ih k n r ) cos sin 2 $, 
+2(s fe + f I {a s cS ) k )b k sin$ r 

-s k ih k b l r sin $ r (sin 2 $» r - |) 
2(s 0fe + f (af ff ),)f fc sin$ 
+ \skih k T l sin$(l + sin 2 $) 



(17) 



(18) 



= q> - \& . 



(15) 



~s u b k T l (2 + 2 cos $ + \ cos $ sin 2 $) (19) 
T r = 2(s Q k + fi(a^ s )k)f k sm^ r 

+ ^s k ih k T l r sm® r {l + sin 2 $ r ) 

~l k ib k T l r (2 + 2 cos $ r + \ cos $ r sin 2 $ r ) (20) 

The usual result from general relativity is obtained for 
vanishing coefficients for Lorentz violation, which follows 
here in the limit B = T = B r = T r = 0. 

The result in Eq. (16) exhibits several interesting fea- 
tures that do not occur in the point-mass limit of GR. 
First, the rotational scalar combinations involving soo, 
Cg , and (affj)o scale the usual GR result. However, due 
to the composition dependence of 5q and (af ff )o, this 
scaling will vary with the central body (or Lens) produc- 
ing the gravitational deflection. Secondly, the anisotropic 
coefficients soj, (sfg)j, and Sjk occur in the deflection 
result such that they are projected along directions asso- 
ciated with the position of the source and reference star 
(n, h r , etc.). This implies that observations made either 
with significantly changing star positions, or made with 
different sets of stars at different locations in the sky, 
will experience different deflections. This effect is inde- 
pendent of the conventional $, <J> r , and dependence 
of the deflection. Thirdly, deflection occurs in the di- 
rections f and r r perpendicular to n and already 



illustrated in Figs. 3. The anisotropic effects imply that 
observations made over time or with many stars would 
yield access to different combinations of coefficients for 
Lorentz violation, thus increasing the potential "parame- 
ter space" of possible types of Lorentz violations to which 
light deflection observations could be sensitive. 



B. Observational analysis 

The form of Eq. (16) indicates that it depends on a 
number of parameters which can be related to specific 
measurable astronomical quantities, in addition to its de- 
pendence on coefficients for Lorentz violation. To illus- 
trate this we describe in this section how analysis might 
proceed. Our post-newtonian coordinate system is taken 
to coincide with the standard Sun-centered celestial equa- 
torial coordinate system adopted for many studies of 
Lorentz violation. The spatial coordinates are denoted by 
X, Y, and Z with X pointing in the direction of the Sun 
at the vernal equinox, and Z is aligned with the Earth's 
rotation axis. The time coordinate is T and is typically 
defined so that T = at the 2000 vernal equinox. De- 
tails of the Sun-centered celestial equatorial coordinate 
system can be found in Refs. [11, 32]. In particular, the 
reader is referred to a depiction of this coordinate system 
in Fig. 1 of Rcf. [11]. 

In general, for data analysis, one can seek to express 
the relative deflection (16) such that its functional form 
is 

5* = 5^(6, 9 r , 0, 4> r , T,s^, c s TT , ...). (21) 

The first four parameters (9,4>) and (9 r ,<fi r ) determine, 
respectively, the direction of h and fi r relative to the 
Sun-centered frame. Explicitly, the unit vectors take the 
form 

h = — (sin 9 cos (j>, sin 9 sin <j>, cos6>), 
n r — — (sm9 r cos <p r , sin# r sin0 r , cos 9 r ). (22) 

Specifically, (9, </>) and (6 r , <f> r ) are taken as co-latitude 
and right ascension on the celestial sphere for the source 
star and reference star, respectively. In terms of these 
angles, the angle *o can be determined using Eq. (14) 
and the angles $ can $ r can be determined using f p -h = 
cos $ and f p ■ h r — cos $ r . 2 

The time parameter T appears in part due to the ob- 
server's motion in the Sun-centered frame. To sufficient 
accuracy in any of the Lorentz- violating terms of (21) 



2 In principle, the coordinates of the source star and reference 
star are also affected by gravitational deflection. However, the 
deflection expression (21) is already expressed at PNO(2), and 
so within this expression, initially measured or "catalog" values 
can be used. 
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we can assume that the motion of the observer is explic- 
itly known. For example, for the Earth we can assume a 
circular orbit and use 

f p = (— cos £1T, — cos rj sin flT, — sin r\ sin flT) , (23) 

where r\ is the inclination of the ecliptic to the equatorial 
plane and fi is the Earth's orbital frequency. 

The last set of parameters s^, (a^ ff ) M , and cf T are the 
coefficients for Lorentz violation expressed in the Sun- 
centered frame coordinates. In addition to the dot prod- 
ucts (n r -b, h r - f , etc.), the terms occurring in the deflec- 
tion formula (16) contain projections of the coefficients 
along the six unit vectors for the source and reference 
source {n, b, f, h r , b r , f r }. To capture the orientation de- 
pendence of the coefficients we must express all of these 
unit vectors in the Sun-centered frame. This can be ac- 
complished using (22) and (23). The impact parameter 
vector is defined by b = r p — n(n ■ r p ) (see Fig. 1). Using 
this, the unit vectors b and b r can be written as 

t f v — n cos $ 

b = ■ rh ' 

sin <P 

b r = f P-*r (24) 

sin$ r 

From this expression and (22) the unit vectors f and f r 
can be constructed using the cross products 

f = n x b, 

r r = h r x b r . (25) 

The full set of unit vectors allows us to express projec- 
tions of the coefficients in terms of Sun-centered frame 
quantities. For example, if only s X z happens to be 
nonzero, we obtain for the projection ~§jkt j I > , 

SjKT J b K = s xz (r x b z +r z b x ) 
= esc 2 $sjcz[sin^ 

x (sin <f) cos OT — cos rj cos cj> sin Q,T) 
x (cos QT — sin 9 cos $ cos <p) 
+ sin f2T(sin 77 sin 9 sin <f> — cos 77 cos 9) 
x (cos 6 cos $ - sin rj sin QT)] , (26) 

where we assumed the case of an Earth observer and 
made use of the formula (23) for f p . 

The full expression of the deflection angle (16), in 
terms of the parameters indicated in (21), can be ob- 
tained using the methods just described. This calculation 
is lengthy, and so we do not include it here, but it should 
be straightforward to include in a suitable data analysis 
code. To give a flavor of the leading Lorentz-violating 
effects in light bending, we determine the approximate 
form of (16) for a special case in the next subsection. 

Finally we note that we have included ellipses in ex- 
pression (21) to indicate that we have not discussed a 
number of other astronomical parameters that may be 
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important to a rigorous analysis, such as the effects of 
parallax and proper motion. Furthermore, we have ne- 
glected aberration effects in the result (16) that depend 
on the Earth's or observer's velocity. We find that one 
effect of these terms is that they multiply the coefficients 
for Lorentz violation in (16) by higher powers of velocity, 
resulting in PNO(3) effects. In addition, terms propor- 
tional to v or its square that arise for aberration in the 
conventional case [33], can implicitly depend on the co- 
efficients for Lorentz violation through orbital dynamics. 
The effects of Lorentz violation on orbits can in principle 
be incorporated using the equations of motion derived in 
Rcfs. [13, 15]. 



C. Conjunction example 

To elucidate the Lorentz-violating effects in the main 
result (16), we work with a special case of the general 
SME deflection result that involves only one set of coeffi- 
cients for Lorentz violation. We set to zero all coefficients 
save those contained in sjk- We also ignore any contri- 
butions from sjj — stt that scale the GR deflection. 

We focus in this example on an Earth observer during 
times near the summer solstice when the Earth lies be- 
low the negative Y axis in the Sun-centered frame (the 
summer solstice occurs when SIT = n/2). To exploit the 
peak behavior of the deflection result (16), we suppose 
that the source star is located on the celestial sphere so 
that its light just grazes the Sun on its way to earth. The 
reference star is taken to be a considerable angular dis- 
tance away from the source star so that its gravitational 
deflection can be ignored to a good approximation. For 
simplicity, both stars are assumed to have <fi = n/2 = <p r . 

In this near-conjunction scenario, the formula for the 
unit vector in the direction of the Earth observer is given 
by (23) while the unit vector for the source star is given 

by 



n = — (0, cos(?7 + e), sin(?7 + e)), 



(27) 



where e is a small angle indicative of the near-conjunction 
approximation adopted. The other needed unit vectors 
b and f can be obtained from the results in the previous 
subsection. 

If we focus on only the dominant terms in this sce- 
nario, the result (16) simplifies considerably. As pre- 
viously stated, the second term in (16) proportional to 
GM/b r that involves the reference star quantities can 
be neglected since b « b r . We can also discard any 
terms with one or more powers of sin<f>, since they will 
be suppressed in this limit. Finally, using the small angle 
approximation for e can further simplify the expression. 

If t is the time measured from the summer solstice, the 
GR portion of the deflection becomes 



5f 



GR 



-2- 



GMr. 



sine 



R (1 — cos Sit cose) 



(28) 



where R is the earth's orbital radius and Mq is the 
Sun's mass. The portion of the deflection controlled by 
the coefficients for Lorentz violation ~sjk stems from the 
sjKT J b term in Eq. (16). It is given by 



(5* 



GMr. 



LV 







sinOi 



R (1 — cos Sit cose) 2 
x (sa sin Sit cos Sit sin e 

+ s B [cos 2 msin 2 e-sin 2 m]), (29) 

where the combinations of coefficients sa and sg are 
given by 



sa 

SB 



sin 2 T)Syy + cos 2 rjszz — sxx — sm2r)SYZ, 



smr]SxY — cosrjsxz- 



(30) 



The result from GR in (28) is plotted in Fig. 5 near the 
time of conjunction (t = 0) using the values r] — 23.4° 
and e = 0.27° (grazing limit). This deflection is peaked 
and is symmetrical around t = 0. For the Sun as the 
deflecting body, the peak value of 1.75" is well known. 
Note that the sign of the GR deflection is negative. This 
is consistent with the (outward) apparent deflection de- 
picted in Fig. 2, since i5\1/gr = * — *o is the difference in 
the observed angle from the unperturbed, or zeroth-order 
angle. 



(arcsec) 
0.5 r 




t (days) 



FIG. 5: The behavior of the deflection Sty near conjunction 
plotted as function of time t around the summer solstice as 
the Sun moves across the field of view. The solid curve is 
the deflection from GR, the dashed curve is the deflection 
amplitude from the coefficients sa, and the dotted curve is 
the deflection amplitude from the coefficients sb- 



The two types of Lorentz-violating signals in (29) are 
also plotted in Fig. 5 near the time of conjunction using 
the same assumptions on i] and e. For these curves we 
plot amplitudes, or "partials" , for each of the coefficient 
combinations sa and sb- It is evident that these sig- 
nals arc qualitatively different from the GR case. The 
amplitude for sa displays symmetrical behavior around 
t = while the amplitude for sb is antisymmetric around 
t = 0. Both of these signals display mildly oscillatory be- 
havior near t = that could potentially be useful for 
analysis. 
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If measurements were obtained for differing orien- 
tations of the observer, the signals for Lorentz viola- 
tion would involve coefficient combinations distinct from 
those in Eq. (30). As an example of this orientation de- 
pendence, suppose instead that the (near-conjunction) 
observations took place with the Earth near the vernal 
equinox, when the Sun is along the positive X axis of the 
chosen coordinate system. The GR result is very similar 
to Eq. (28) in this case, but the Lorentz-violating piece 
(5^lv differs in its details. Specifically, for this configu- 
ration we find 

^ _ GMq cos 7] sin ilT 

LV R (1 — cos QT cose — sin r/ sine sin SIT) 2 

x [{syy - szz){ \ sin2?7sin 2 SIT 

— cos i] sin SIT cos SIT sin e) 
+s Y z (cos 2 SIT sin 2 e - cos 2i] sin 2 SIT 

- 2 sin r) sin e sin SIT cos SIT)] , (31) 

where T is measured from the conjunction time T = 0. 
It is evident from this result that the signal depends on 
the coefficients combinations syy — szz and syz, rather 
than those given in Eqs. (30). 

IV. SUMMARY AND ESTIMATES 

In this work we have identified the dominant signals for 
local Lorentz violation in light-bending observations. A 
general formula making use of euclidean arclength that 
is valid for the post-newtonian limit of any stationary 
metric was established in (5). Working within the SME 
effective field theory framework, we applied this formula 
to the deflection of a light ray from a straight line path in 
the grazing limit in Eq. (9), and more accurately in Eq. 
(10). The results display anisotropic behavior of light 
bending controlled by coefficients for Lorentz violation, 
as demonstrated in Fig. 3. 

In the latter part of this work, we calculate a more 
practical formula for the change in the measured angle 
* between two stars in Eq. (16). We describe generally 
how this result can be expressed in terms of astronomical 
quantities suitable for data analysis. The approximate 
behavior of the deflection angle shift 5^ on the coeffi- 
cients for Lorentz violation was elucidated with a spe- 
cific near-conjunction example. We compare this uncon- 
ventional behavior with the standard behavior predicted 
from GR (see Fig. 5). 

It would be of interest to perform a rigorous analysis of 
potential sensitivities for future missions, perhaps involv- 
ing detailed simulations. Such simulations have already 
been performed for the GR and PPN case for the planned 
Gaia mission [31, 35]. The starting point for such sim- 
ulations is typically the relative deflection between two 
stars. We provide this expression for the SME in Eqs. 
(12) and (16) of this work. A numerical least squares es- 
timation could be attempted using the partial derivatives 
of <3/ with respect to the coefficients for Lorentz violation 



s^v, (afffV, an d Cq , along with other relevant parame- 
ters. 



Observatory 




s T j + {a/M){a s cS )j 


SJK 


Ref. 


LATOR 


1(T 8 


1(T 8 


lO" 7 


[37] 


Gaia 


icr 6 


icr 6 


1(T 5 


[31] 


Hipparcos 


lO" 3 


1(T 3 


icr 2 


[36] 


Optical 


HT 1 


1 


i 


[3] 



TABLE I: Crude estimates of sensitivities of current and 
future light-bending tests to various combinations of coeffi- 
cients for Lorentz violation. The shorthand defined by 
s' TT = stt + c TT + If (of ff ) T . 

Though it is beyond the scope of this work to per- 
form detailed simulations for future missions or analysis 
of available data from current and past observations, we 
provide in Table I order of magnitude estimates of sen- 
sitivities to different combinations of coefficients. These 
estimates are based on existing constraints on deviations 
from GR in light-bending tests or projected measurement 
accuracies. We include the proposed Laser Astromet- 
ric Test of Relativity (LATOR) [37], the planned Gaia 
mission [34], the past Hipparcos mission [36], and past 
ground-based optical observations [3] . This is not a com- 
prehensive list and other dedicated light-bending tests 
may also be of interest [38]. 

The estimates in Table I can be contrasted with ex- 
isting constraints on coefficients for Lorentz violation. 
The rotational scalar combination stt has not been for- 
mally constrained by rigorous data analysis. Care must 
be taken in light-bending tests since this coefficients also 
appears in the Newtonian force law multiplying the com- 
bination GM . Therefore stt is expected to be correlated 
with orbital tests. Combining light-bending results with 
results from orbital tests could yield the measurements of 
stt at the levels indicated in the table. Similar consider- 
ations hold for the scalar matter coefficients (afff)r and 
c TT . Details on this type of comparative measurement 
can be found in Refs. [14, 15]. Note that the matter coef- 
ficients can in principle be separated from stt by using 
deflecting bodies of differing composition. Alternatively, 
one can combine results from light-bending observations 
with current constraints on the matter sector coefficients 
from earth laboratory experiments [19]. 

The coefficient combination stj + (a/M)(af ff )j is also 
of primary interest for light-bending tests. The three co- 
efficients stj are currently constrained at the 10~ 5 — 10~ 6 
level from recent lunar laser ranging and atom interfer- 
ometry tests [11]. For a source body composed of or- 
dinary matter, the source combination a(a^ s )j depends 
on the coefficients for the electron e, proton p, and neu- 
tron n. Specifically for the Sun as the source body, we 
have a{a s eS )j/M « 0.5GeV- 1 a[(a| ff ) J + (^ ff ) J + (aS r ) J ] 
[15]. Future missions, such as Gaia or LATOR, could 
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tighten the constraints on stj and perform the first anal- 
ysis of astrophysical measurements of the matter sector 
coefficients {a% s )j, {a^ s )j, and (a™ s )j at the competitive 
1CT 6 GeV level or better. 

Other related tests are also of interest. This includes 
tests involving the classic time-delay effect [39] . As men- 
tioned previously in this work, modifications to the time 
delay formula arise from local Lorentz violation and have 
been analyzed in Refs. [13, 15]. Note that the signals for 
Lorentz violation in the time-delay effect and light bend- 
ing differ in their dependence on coefficients for Lorentz 



violation, as discussed in Sec. II C. It is therefore of in- 
terest to consider all such tests, as they could be used to 
place independent constraints on Lorentz violation. 

Finally we note that we have not treated here the broad 
subject of gravitational lensing. The deflection angle for- 
mulas calculated in this work could form a starting point 
for analysis [17]. A comprehensive investigation of the ef- 
fects of local Lorentz violation on weak and strong grav- 
itational lensing would be of definite interest but lies be- 
yond the scope of this work [40] . 
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